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Abstract 

Background: Physician collaboration, which evolves among physicians during the course of providing healthcare 
services to hospitalised patients, has been seen crucial to effective patient outcomes in healthcare organisations 
and hospitals. This study aims to explore physician collaborations using measures of social network analysis (SNA) 
and exponential random graph (ERG) model. 

Methods: Based on the underlying assumption that collaborations evolve among physicians when they visit a 
common hospitalised patient, this study first proposes an approach to map collaboration network among 
physicians from the details of their visits to patients. This paper terms this network as physician collaboration 
network (PCN). Second, SNA measures of degree centralisation, betweenness centralisation and density are used to 
examine the impact of SNA measures on hospitalisation cost and readmission rate. As a control variable, the impact 
of patient age on the relation between network measures (i.e. degree centralisation, betweenness centralisation and 
density) and hospital outcome variables (i.e. hospitalisation cost and readmission rate) are also explored. Finally, ERG 
models are developed to identify micro-level structural properties of (i) high-cost versus low-cost PCN; and 
(ii) high-readmission rate versus low-readmission rate PCN. An electronic health insurance claim dataset of a very 
large Australian health insurance organisation is utilised to construct and explore PCN in this study. 

Results: It is revealed that the density of PCN is positively correlated with hospitalisation cost and readmission rate. 
In contrast, betweenness centralisation is found negatively correlated with hospitalisation cost and readmission rate. 
Degree centralisation shows a negative correlation with readmission rate, but does not show any correlation with 
hospitalisation cost. Patient age does not have any impact for the relation of SNA measures with hospitalisation cost 
and hospital readmission rate. The 2-star parameter of ERG model has significant impact on hospitalisation cost. 
Furthermore, it is found that alternative-k-star and alternative-k-two-path parameters of ERG model have impact on 
readmission rate. 

Conclusions: Collaboration structures among physicians affect hospitalisation cost and hospital readmission rate. The 
implications of the findings of this study in terms of their potentiality in developing guidelines to improve the 
performance of collaborative environments among healthcare professionals within healthcare organisations are 
discussed in this paper. 

Keywords: Physician collaboration network, Exponential random graph, Social network analysis, Hospitalisation cost 
and readmission rate 
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Background 

Collaborations among physicians have been found very 
important to the effectiveness in delivering healthcare 
services and in producing better patient outcomes [1,2]. 
The structure of collaboration (i.e. the way how people 
communicate and collaborate with others in a collabora- 
tive environment) among the hospital staff could not be 
the same in different hospitals or healthcare organisations. 
These various structures may have different impact on 
healthcare outcome measures (e.g. hospitalisation expenses 
and patient satisfaction) in various healthcare contexts [3]. 
Some structures could be more conducive in terms of pa- 
tient and hospital outcomes compared to others. Therefore, 
it is necessary to analyse different structures of collabora- 
tions among healthcare professionals and their impact on 
outcome variables. In this paper, measures of social net- 
work analysis (SNA) and exponential random graph (ERG) 
models are employed to explore physician collaborations in 
order to find out structural attributes of physician colla- 
borations that are conducive to hospitalisation cost and 
readmission rate. 

Collaboration, which is a recurring process where two 
or more people or organisations work together towards 
common goals [4], enables individuals and organisations 
to work together more effectively and efficiently. Colla- 
borative relationships among individuals are highly cele- 
brated in organisations because the synergies realised by 
combining multi-dimensional efforts and diverse expert- 
ise produce benefits greater than those achieved through 
individual effort [5]. In the context of healthcare service 
providers or hospitals, collaboration among different 
healthcare professionals is recognised as a catalyst to 
improved patient outcomes such as less hospital length 
of stay and hospitalisation cost [6-8], lower death rate 
[9] and higher satisfaction [10,11]. In healthcare settings, 
collaboration allows input from multiple professions (e.g. 
nurse and physicians), which could produce decisions 
leading to better patient outcomes because those decisions 
are based on more complete information [8]. 

The context of this study is the physicians' collabo- 
rations that evolve within healthcare service providers or 
hospitals during the course of providing healthcare ser- 
vices to patients. Arguably, it can be conceptualised that 
physicians collaborate with each other and with other 
hospital staffs (e.g. nurses) in order to provide effective 
services to hospitalised patients. Based on the patient 
condition and unavailability of their colleagues, phy- 
sicians might seek advices or suggestions from other 
physicians working in different workplaces. Because of 
this type of medical practice culture in healthcare service 
providers or hospitals, a professional collaboration net- 
work has eventually been developed over time among 
physicians. This study terms this network as 'Physician 
Collaboration Network (PCN)'. 



The measures and methods of social network analysis 
(SNA) have been found useful in investigating networks 
(e.g. PCN) and their effects on performance [12,13]. 
SNA can be seen as the mapping and measuring of 
relationships among participating actors [14] and can 
provide both a visual and a mathematical analysis of 
network relations among actors. It plays an important 
role in identifying and quantifying the informal network 
which functions at a level beyond the formal and trad- 
itional organisational structure of actor relationships 
[13]. In modelling structures of PCN, this study uses the 
Exponential Random Graph (ERG) model which is a 
probabilistic model and has been utilised extensively in 
the social science literature to study the dynamics of 
network formation from underlying locally prominent 
micro structures such as 2-star, 3-star, triangle and so on 
[15]. Although most of the studies about ERG focus on 
building the theory of ERG models, recently researchers 
have applied ERG models in practice, such as, to under- 
stand whether external connections beyond the department 
are important to the understanding of the departmental 
structure of an Australian Government Organisation [16], 
to explore the dynamics of biological networks [17] and to 
examining the communication dynamics of networks 
under stress [18]. This study considers hospitalisation cost 
and readmission rate as surrogate measures for the effect- 
iveness and efficiency of physician collaborations. There 
are several evidences of the use of readmission rate and 
hospitalisation cost as outcome measures in the healthcare 
literature [19-21]. Further, this study utilises patient age as 
a control variable. The use of patient age as a control or 
moderating variable has been found in several studies of 
the present healthcare literature [22,23] . 

The aim of this study is threefold. It first proposes a way 
to map physicians' collaboration from their visiting infor- 
mation to patients. Then this study explores, by considering 
patient age as control variable, what macro-level (i.e. the 
complete structure of a collaboration network) SNA mea- 
sures of these collaboration networks affect hospitalisation 
cost and hospital readmission rate. The last aim of this 
study is to examine what micro-level structures (i.e. small 
structures among few physicians in a PCN such as 2-star) 
among physicians affect hospitalisation cost and hospital 
readmission rate. For this purpose, this study considers the 
only top 5 collaboration networks having low-cost and low- 
readmission rate, and top 5 collaboration networks having 
high-cost and high-readmission rate. The organisation of 
this paper is as follows. The rest of the 'background' section 
reviews the current collaboration literature in healthcare 
context and illustrates the way to map physicians' collabor- 
ation from their visiting information to patients. The re- 
search methodologies (i.e. description of SNA measures, 
ERG models, research dataset, dependent variables and 
control variable) followed in this study are described in the 
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'methods' section. The 'results' section describes the findings 
of this study. Finally, the 'discussion and conclusion' section 
discusses the findings of this study and provides some policy 
recommendations for healthcare managers or administra- 
tors. This section also makes a conclusion for this paper. 

Literature review: collaboration in healthcare context 

There are numerous studies in current literature exploring 
the effect of collaboration among healthcare professionals 
on patient outcomes and hospital performance. Most of 
these studies explore hospital performance and patient 
outcomes by analysing collaboration networks among 
different healthcare professionals such as nurse-physician 
collaboration [9], physician-pharmacist collaboration [24], 
physician-patient collaboration [25], hospital-physician 
collaboration [26], and inter-professional and interdis- 
ciplinary collaboration [27]. Cunningham et al. [28] have 
conducted an orderly review of studies of professionals' 
network structures, analysing factors connected with 
network effectiveness and sustainability, specifically in re- 
lation to the quality of care and patient safety. The authors 
explore MEDLINE, CINAHL, EMBASE, Web of Science 
and Business Source Premier from January 1995 to 
December 2009. A majority of 26 studies reviewed used 
social network analysis to analyse structural relationships 
in networks: structural relationships within and between 
networks, health professionals and their social context, 
health collaborations and partnerships, and knowledge 
sharing networks. Essential features of networks explored 
were administrative and clinical exchanges, network per- 
formance, integration, stability and influence on the qual- 
ity of healthcare. They have also noticed that more recent 
studies demonstrate that cohesive and collaborative health 
professional networks can promote the coordination of 
care and contribute to improving quality and safety of 
care. Structural network vulnerabilities include cliques, 
professional and gender homophily and over-reliance on 
central agencies or individuals. Efficient professional 
networks engage basic structural network features (e.g. 
bridge, broker, density, centrality, degree of separation, 
social capital and trust) in generating collaboratively 
oriented healthcare. This requires effective transmission of 
information and social and professional interaction within 
and across networks. For those using networks to improve 
care, recurring success factors are: understanding your 
networks characteristics, attending to its functioning and 
investing time in facilitating its improvement. Despite this, 
there is no guarantee that time spent on networks will 
necessarily improve patient care. 

Another classic study, led by Knaus and his team, identi- 
fies a significant relationship between the degree of nurse- 
physician collaboration and patient mortality in intensive 
care units [9]. They study treatment and outcome in 5030 
intensive care unit patients and find that hospitals where 



nurse-physician collaboration is presented report a lower 
mortality rate compared to the predicted number of pa- 
tient deaths. Conversely, hospitals that are noted for poor 
communication among healthcare professionals exceed 
their predicted number of patient deaths. In a two group 
quasi-experiment on 1207 general medicine patients (n = 
581 in the experimental group who received care from a 
specially designed care management plan that facilitated 
higher collaboration among hospital staff and n = 626 in 
the control group who received the usual care), Cowan 
et al. [6] notice average hospital length of stay, total hospi- 
talisation cost, and hospital readmission rate are signi- 
ficantly lower for patients in the experimental group than 
the control group (5 versus 6 days, P < .0001) which 
contributes a ( bacl<fill profit' of US$1591 per patient to 
hospitals. There are other studies that also highlight the 
importance of collaboration among healthcare profes- 
sionals for better patient outcomes. 

Sommers et al. [29] examine the impact of an interdis- 
ciplinary and collaborative practice intervention involving 
a primary care physician, a nurse, and a social worker for 
community-dwelling seniors with chronic illnesses. They 
conduct a controlled cohort study of 543 patients in 18 
private office practices of primary care physicians. The 
intervention group receives care from their primary care 
physician working with a registered nurse and a social 
worker, while the control group receives care as usual 
from primary care physicians. They notice that the inter- 
vention group produced better result to readmission rate 
and average office visits to all physicians. Moreover, the 
patients in the intervention group report an increase in 
social activities compared with the control groups de- 
crease. There are other studies emphasising collaboration 
for effective patient outcome across professional bounda- 
ries within hospitals. By analysing data collected from 105 
interviews (with 40 physician, 32 case managers, 23 phy- 
sician office staff, 8 administrators, and 2 case assistants), 
Netting and Williams [30] argue that there is a growing 
need to collaborate and communicate across professional 
lines rather than make assumptions about who can do 
what for better patient outcomes, professional satisfaction, 
and hospital performance. 

Like these studies, most of the collaboration studies of 
contemporary healthcare literature advocate for the effec- 
tive and efficient collaboration among healthcare pro- 
fessionals for better patient care. Proper collaborations 
among hospital staff positively drive both total hospitalisa- 
tion cost and hospital readmission rate. There are many 
other studies in healthcare context that analyse networked 
collaboration among healthcare specialists to explore 
different aspects of professional behaviour and quality 
patient care, such as, to evaluate the effects of GP network 
organisation on their prescribing behaviour [31] and to 
develop a selection criteria of group members in order to 
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improve the effectiveness of team-based approach to 
patient care [32]. However, none of these studies, to our 
knowledge, provides any guidelines: (i) about the network 
structure of effective collaboration; (ii) what type of collab- 
oration structure is more conducive compared to others; 
and (iii) on how individual healthcare professional should 
develop relations with others over time in a collaborative 
environment for better performance. This study considers 
only the physician collaboration and addresses all of these 
three issues by exploring PCNs using measures of social 
network analysis (SNA) and exponential random graph 
(ERG) models. 

Physician collaboration network (PCN) 

Collaboration in healthcare is defined as healthcare 
professionals accepting complementary roles and jointly 
working together, sharing accountability for problem- 
solving and making decisions to develop and implement 
plans for patient care [33,34]. Collaboration among 
physicians, nurses and other healthcare professionals 
increases team members' perception of each other s type 
of knowledge and skills, leading to continued improve- 
ment in decision-making [35]. It can take place in both 
face-to-face interactions and electronically via fast-paced 
encounters such as e-mail. In whatever location or form, 
collaboration includes an exchange of beliefs and ideas 
that acknowledges the perspectives of all collaborators, 
whether or not agreement is accomplished in the inter- 
action [36]. To minimise misunderstandings, it is also 
essential to define what is not implied by the term 
'collaboration . It does not imply supervision, nor is it 
simply a one-way or two-way information exchange. 
Efficient professional collaborative relationships require 
mutual respect [37]. They also call for trust and dili- 
gence. In complex and sophisticated healthcare systems, 
collaboration is generally challenging. Collaboration 
may seem idealistic and perhaps even non-realistic. 
However, Kramer and Schmalenberg [37] state that 
collaborative partnerships are worth the effort because 
they result in improved effects for patients as well as 
individual development for collaborators. 

Collaboration between physicians has been poorly inves- 
tigated; the overwhelming focus of research on physicians 
has been on their collaboration with patients [38]. Inter- 
estingly, research on collaboration between physicians has 
focused on the discussion of medical mistakes, collegial 
control and other negatively recognised aspects of medical 
care [39-42]. One area that has been seen critical is the 
culture of medicine and the socialisation of medical 
students, interns, and residents into that culture by physi- 
cians [43,44]. Atkinson [38] argues that "biomedical know- 
ledge is socially produced and culturally specific . . . [and] 
dependent upon certain fundamental features of medical 
culture, which is itself produced and reproduced through 



processes of socialisation (p. 46). His study of haema- 
tologists' consultations with physicians of other specialties 
demonstrates the method of generating medical know- 
ledge through collaboration among physicians. According 
to Atkinson, physicians' communication is not the way to 
the accomplishment of medical work; the communication 
is the work. 

In this study, it is assumed that collaborations among 
physicians emerge when they visit common hospitalised 
patients. It is a standard professional practice around the 
world that when physicians visit patients, they give advice 
or suggestions to patients based on their health condition 
and previous medication history deposited in the patient 
log book. All previous advice or suggestions by any 
physician to a patient have been taken into consideration 
during any subsequent physician visit to that patient. This 
kind of practice culture in healthcare organisations or hos- 
pitals enables us to map and, eventually, to model PCNs. 

When physicians visit common patients within the same 
hospital or healthcare organisation PCN emerges among 
them. Figure 1 illustrates an example of such a PCN con- 
struction. In a hospital (say HI), patient Pal is visited by 
Phi, Ph2 and Ph4 physicians, and patient Pa2 is visited by 
Ph2, Ph3 and Ph4 physicians, and physician Ph3 and Ph4 
visit patient Pa3. This is depicted in the patient-physician 
network in Figure 1(a). The corresponding PCN for this 
patient-physician network is demonstrated in Figure 1(b). 
In this PCN, there are network connections with weight 1 
between Phi and Ph2, between Phi and Ph4, and between 
Ph2 and Ph3 because they visit only one common patient. 
The weight of the links between Ph2 and Ph4, and between 
Ph3 and Ph4 are 2 as they have two common patients. 

Methods 

In explaining different SNA measures and ERG model, 
this study uses the terms actor(s) and node(s) interchange- 
ably. Similarly, the words link(s) and tie(s) are exchange- 
able in this paper. 

Measures of social network analysis (SNA) 

This study utilises SNA measures of degree centralisa- 
tion, betweenness centralisation, and network density. 
The selection of these three measures is guided by 
two network theories: (i) Bavelas' Centralisation Theory 
[45]; and (ii) Freeman's Centrality Theory [46]. These 
two theories can explain structural influences of colla- 
boration and communication networks on the group 
performance. 

Degree centralisation and betweenness centralisation 

Before explaining degree centralisation and betweenness 
centralisation, it is required to define degree centrality 
and betweenness centrality. Centralisation is a network- 
level measure whereas centrality is a node-level measure; 
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Patient-physician network 



(b) Corresponding PCN 



Figure 1 Conceptualisation of the collaboration network among physicians, (a) Patient-physician network, and (b) Corresponding PCN (Pa 
stands for patient, and Ph stands for physician. 



thus, the later one needs to be explained first before 
describing the former one. Centrality is an important con- 
cept in studying networks. In conceptual terms, centrality 
measures how central an individual is positioned in a 
network. Degree centrality is one of basic measures of net- 
work centrality. For an actor, it is the proportion of nodes 
that are adjacent to that actor in a network. It highlights 
the node with the most links to other actors in a network, 
and can be defined by the following equation for the actor 
i in a network having N actors [13]: 



d(nj) 
N-l 



(1) 



Where, the subscript D for degree and d(n^ indicates 
the number of actors with whom actor i is connected. 
The maximum value for C D (ni) is 1 when actor i is 
linked with all other actors in the network. For an isolate 
actor, its value is 0. 

Betweenness centrality views an actor as being in a 
favoured position to the extent that the actor falls on the 
shortest paths between other pairs of actors in the net- 
work. That is, actors that occur on many shortest paths 
between the other pair of nodes have higher betweenness 
centrality than those they do not [46]. The betweenness 
centrality for an actor n t (i.e. C B (ni)) can be represented 
by the following equation [13]: 



E 



£/*("<) 



[(AT-l)(JV-2)]/2 



(2) 



Where, i * j * k; gjkind represents the number of the 
shortest paths linking the two actors that contain actor i; 
and gjk is the number of the shortest paths linking actor 
j and k. For the central actor of a start, d B (rii) will take 
its highest value of 1; however, for any peripheral actor 
of a star C B {n\) will take its minimum value of 0. 

A centralisation measure quantifies the range or vari- 
ability of individual actor indices. The set of degree 
centralities, which represents the collection of degree indi- 
ces of N actors in a network, can be summarised by the 



following equation to measure network degree centralisa- 
tion [47]: 



C D 



J2[C D (n*yC D (nt)] 



(3) 



Where, {C D (nj)} are the degree indices of N actors and 
C D (n*) is the largest observed value in the degree indices. 
For a network, degree centralisation (i.e. the index C D ) 
reaches its maximum value of 1 when one actor chooses 
all other (N-l) actors and the other actors interact only 
with this one (i.e. the situation in a star graph). This 
index (i.e. C D ) attains its minimum value of 0 when all 
degrees are equal (i.e. the situation in a circle graph). 
Thus, C D indicates varying amounts of centralisation of 
degree compared to both star and circle graph. 

Similarly, the set of betweenness centralities, which rep- 
resents the collection of betweenness indices of N actors in 
a network, can be summarised by the following equation 
to measure network betweenness centralisation [47] : 



C B 



j^[c B {n*)-C B { ni ) 

i—l 

{N-l) 



(4) 



Where, {Cg(^)} are the betweenness indices of N actors 

and C B (n*) is the largest observed value in the betweenness 
indices. Freeman [46] demonstrates that betweenness cen- 
tralisation reaches its maximum value of 1 for the star 
graph. Its minimum value of 0 occurs when all actors have 
exactly the same betweenness index. 

Network density 

The density of a network represents the proportion of 
existing ties (or, links) relative to the maximum number of 
possible ties among all actors of that network [13]. The 
density value for a network is 1 only when all the actors of 
that network are connected with each other. On the other 
hand, for a completely sparse network, the density value is 
0, which indicates there is no link exists between any two 
actors of that network. For an undirected network of size 
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N (i.e. have N actors), theoretically there are [N* (N- 1)]/ 
2 (i.e. N C 2 ) possible links among its N actors. If there are 
N t links among its N actors in that network, then, math- 
ematically, density can be defined as [13]: 

2*N t 

Density= ^i) (5) 

Exponential random graph (ERG) models 

ERG model can effectively identify structural properties in 
social networks [48]. This theory-driven modelling ap- 
proach also allows to test the significance of structural 
parameters in the process of the formation of a given 
network [18,49]. For instance, a given cost effective PCN 
may be explored using ERG model to examine what micro 
structures play a statistically significant role in the devel- 
opment process of this PCN. It simplifies a complex struc- 
ture down to a combination of basic parameters. The 
advantage of this approach is that it is very general and 
scalable as the architecture of the graph is represented by 
locally determined explanatory variables, and the choice 
of explanatory variables is quite flexible and can be easily 
revised. The disadvantage of this approach is the difficulty 
in estimating the execution time. The reason for that is 
that ERG models are based on simulation and execution 
time for simulation is always unpredictable [15]. Another 
disadvantage of ERG models is the complex interpretation 
when multiple parameters are considered and the diffi- 
culty to get convergence sometimes [15]. 

This paper follows the notation and terminology 
described in Robins et al. [50]. For each pair i and / of a 
set of N actors, X t j is a network tie variable with X t j = 1 
if there is a network tie from i to and X t j = 0 otherwise. 
This paper specifies %ij as the observed value of Xy with 
X the matrix of all variables and x the matrix of ob- 
served ties of the network. X may be directed or non- 
directed. A configuration is a set of nodes and a subset 
of ties among them. For example, an edge is a subset of 
two nodes in which one node is connected by a tie to 
other, and a 3-star is a subset of four nodes in which 
one node is connected by a tie to each of the other three 
nodes. Similarly, n-star is a subset of n nodes in which 
one node is connected by a tie to each of the other (n-1) 
nodes. Configurations are defined hierarchically, so that a 
triangle also includes three 2-stars. The general form of 
the class of (homogeneous) ERG models is as follows [50] : 

Pr(X = x) = iexp{£ A?Ag »} (6) 

Where, (i) the summation is over configuration types 
A; different sets of configuration types represent differ- 
ent models (e.g. dyadic independence or Markov random 
graph); (ii) rj A is the parameter corresponding to a 



configuration of type A; (iii) g A (x) is the network statis- 
tic corresponding to configuration A (for homogeneous 
Markov graph models this is the number of configura- 
tions of type A observed in the network: for example, 
the number of triangles); and (iv) k is a normalising 
quantity to ensure that Eq. (6) is a proper probability 
distribution. 

A commonly used sub-class of ERG models is the 
Markov random graph in which a possible tie from i to ; 
is assumed conditionally dependent only on other possible 
ties involving i and/or / [51]. This sub-class of ERG model 
is also known as the low-order model. An example of a 
Markov random graph model for non-directed networks, 
with edge (or, density), 2-stan 3-star and triangle parame- 
ters, is given below [52]: 

Pr(X = x) = - exp {6L(x) + a 2 S 2 (x) + a 3 S 3 (x) + tT(x)} 

(?) 

In Eq. (7), 6 is the density or edge parameter and L(x) 
refers to the number of edges in the graph x; o k and 
Sk(x) refer to the parameter associated with k-star effects 
and the number of k-stars in x; while r and T(x) refer to 
the parameter for triangles and the number of triangles, 
respectively. For a given observed network x, parameter 
estimates indicate the strength of effects in the data. For 
instance, a large and positive estimate for o* 2 suggests 
that, given the observed number of edges and stars, net- 
works with more 2-stars are more likely. The configura- 
tions and parameters of Markov random graph model 
(i.e. low-order model) is shown on Figure 2a. These pa- 
rameters relate to some well-known structural regularity 
in the network literature and represent structural ten- 
dencies in the network (e.g. mutuality and transitivity). 
They were chosen because they are conceptualised as 
forces which drive the formation of the network itself. 
For example, transitivity is conceptualised as a force 
which drives the formation of the network itself (the 
friends of our friends are more likely to be our friends). 
Snijders et al. [15] later propose three new configura- 
tions (i.e. alternating k-stars, alternating k-triangles and 
alternating independent two-paths) that can be included 
in specifications for ERG models. They define a new sub- 
class of Markov random graph model (i.e. high-order 
model), which considers parameters of both Figure 2a,b. 
This study utilises both high- and low-order ERG models 
for modelling PCNs. 

There are two methods commonly used in the statistics 
and social network communities to estimate the max- 
imum likelihood fit to ERG models: Markov chain Monte 
Carlo maximum likelihood estimation and maximum 
pseudo-likelihood estimation. They can also be used for 
network simulation. To date, the most common form of 
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estimation for Markov random graph models is the 
maximum pseudo likelihood [53]. The properties of the 
pseudo-likelihood estimator are not well understood and 
the pseudo-likelihood estimates can at best be thought of 
as approximate. Also, it is not clear from existing research 
as to when pseudo-likelihood estimates may be acceptable. 
Therefore, Monte Carlo Markov chain maximum 
(MCMC) likelihood estimation, when available, is the 
preferred estimation procedure. This study also uses this 
estimator. That means this study utilises Markov random 
sub-class of ERG models for modelling cross-sectional 
PCNs and maximum Monte Carlo Markov chain max- 
imum (MCMC) estimator for estimation purpose. 

Research dataset 

This research utilises health insurance claim dataset to 
explore physician collaborations using measures of SNA 
and ERG models. This dataset is provided by a non- 
profit health insurance organisation (i.e. Hospital Con- 
tribution Fund, HCF), which is the third largest health 
insurance organisation in Australia. It includes mem- 
bers' claim data from January 2005 to February 2009. 
This dataset contains mainly three different categories 
of claim information: (i) ancillary claim (lodged by 
hospital); (ii) medical claim (lodged by doctor or physi- 
cians); and (iii) hospital claim (lodged by hospital). 
Ancillary claims are auxiliary claims for medical ser- 
vices such as dental, optical, physiotherapy, dietician, 
and pharmaceutical. All claims lodged by specialist phy- 
sicians, except of the ancillary type, are medical claims. 
The claims for the services provided to hospitalised 
patients in private or public hospitals that are approved 
by the Department of Health, Australia are considered 
as hospital claims. In general, patients have medical 



claims, hospital claims, and very few ancillary claims 
for their admissions to hospitals. This study uses the claim 
information to construct PCN for a particular type of 
hospitalised patients (e.g. knee surgery patient). 

As people have hospital admissions for a wide range of 
illness and patients with a particular disease need to be 
seen by particular specialist physicians, different types of 
PCNs (e.g. a PCN for knee surgery patients and a PCN 
for heart surgery patients) are being evolved inside a 
hospital for hospitalised patients suffering from different 
types of diseases. For research analysis purpose, this 
study considers PCNs only for total hip replacement 
(THR) patients from 85 different hospitals where at least 
5 THR patients get admitted during the data collection 
period. So, 85 PCNs evolved within these 85 hospitals. 
In these hospitals, 2229 patients get admitted during our 
data collection period. These patients lodged in total 
1383 ancillary claims, 65871 medical claims, and 23369 
hospital claims. The basic statistics of these 85 PCNs is 
given in Table 1 (last column). 

To explore physician collaborations using SNA mea- 
sures, these 85 PCNs are used in this research. However, 
for ERG modelling, this study utilises 20 PCNs. In particu- 
lar, this study considers 5 most expensive PCNs (termed 
as high-cost PCNs) and 5 least expensive PCNs (termed as 
low-cost PCNs) to explore, using ERG models, how micro- 
level network structures varied for PCNs having different 
total hospitalisation cost. For modelling PCN in terms of 
readmission rate, this study considers 5 PCNs that have 
the highest readmission rate (termed as high-readmission 
PCNs) and compare their structures with 5 PCNs which 
have the lowest readmission rate (termed as low-readmission 
PCNs). The basic statistics of these 20 PCNs is given in 
Table 1 (the first 5 columns). 




Parameters for low-order models (^0 Additional parameters for high-order models 

Figure 2 Configurations and parameters for exponential random graph models [50]. 

V J 
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Dependent variable: Hospitalisation cost and readmission 
rate 

Hospitalisation cost 

In calculating hospitalisation cost, this study considers all 
payments made by the health insurance organisation for 
each THR patient to the health service providers, regard- 
less of how much that patient pays in return to the health 
insurance organisation (which depends on the health in- 
surance policy type and percentage of coverage amount). 

Readmission rate 

For any PCN, readmission rate represents the ratio of 
patients (in percentage) who have hospital admissions 
more than once for their THR surgeries. That means, 

Readmission Rate (PCN) 

No. of THR patient(s) who admitted more than once 

£ — r^— ; X 100% 

Total admitted THR patients 

(8) 

Control variable: patient age 

For each PCN, the average age for all patients is calcu- 
lated. This average age is considered as control variable 
to investigate whether the patient age has any impact for 
the relation of SNA measures of PCN with hospitalisa- 
tion cost and hospital readmission rate. 

Results 

This section reports the results of this study. 

Mapping physician collaboration network (PCN) from 
insurance claim dataset 

From the medical claim details of HCF dataset, the num- 
ber of physicians visit a particular hospitalised patient 
during her or his hospitalisation period can be revealed 
because physicians make a medical claim to HCF for 
every single visit to hospitalised patients. Based on this 
information and by applying the PCN development ap- 
proach (as illustrated in Figure 1) and process (described 
in section 3), the structure of PCN of each hospital for 
THR patients has been constructed. An example of the 



construction of PCN structure from the research dataset 
is given in Figure 3. Organisation Risk Analyser (ORA), 
which is a meta-network assessment and analysis tool 
[54], is utilised to construct such PCNs. 

Social network analysis (SNA) measure and physician 
collaboration network (PCN) 

SNA measures and their impact on the hospitalisation 
cost and readmission rate are illustrated in Table 2. This 
table also shows the descriptive statistics (i.e. mean and 
standard deviation) of each variable. Although degree 
centralisation does not show a correlation with hospita- 
lisation cost (rho = 0.112, p>0.05 at 2-tailed), it shows 
negative correlation with readmission rate (rho = - 0.373, 
p<0.01 at 2-tailed). An increase in degree centralisation 
produces a downturn for readmission rate. Density of 
PCN is positively correlated with both hospitalisation cost 
(rho = 0.282, p<0.01 at 2-tailed) and readmission rate 
(rho = 0.358, p<0.01 at 2-tailed). Both hospitalisation cost 
and readmission rate of a hospital change proportionally 
with the change in the density of the PCN of that hospital. 
On the other hand, the correlation coefficient values of 
Table 2 reveals that betweenness centralisation of PCN 
is negatively correlated with both hospitalisation cost 
(rho = -0.264, p<0.05 at 2-tailed) and readmission rate 
(rho = -0.283, p<0.01 at 2-tailed). As it is always 
expected to have low hospitalisation cost and readmis- 
sion rate, this result indicates that low betweenness 
centralisation is not conducive for healthcare service 
providers or hospitals. We also develop simple linear re- 
gression models for each of hospital outcome variables 
(i.e. hospitalisation cost and readmission rate) and PCN 
estimates (i.e. degree centralisation, betweenness centralisa- 
tion and density). These models, as described in Table 3, 
allow checking relative influence and independence in 
the associations of independent network variables and 
dependent hospital outcome variables. All models, except 
the first model (i.e. considering degree centralisation and 
hospitalisation cost), show statistically significant output. 

The effect of patient age as control (or moderating) 
variable is summarised in Table 4. We develop regression 
models by considering each of the network measures and 



Table 1 Summary statistics of 5 low-cost, 5 high-cost, 5 low-readmission, 5-high readmission rate and the total 85 PCNs 

Item Hospitalisation cost Readmission rate All PCNs (85) 



Low (5) High (5) Low (5) High (5) 



Average cost per patient ($AUD) 


16582 


29949 


18931 


26400 


24010 


Average readmission rate (%) 


7.47 


10.98 


0.00 


21.94 


11.64 


Average number of patient per PCN 


61.8 


85.8 


19.2 


30.8 


26.22 


Average number of doctors' visit per patient 


17.01 


34.79 


29.55 


26.66 


26.03 


Average patients' age (year) 


54.08 


67.87 


73.06 


65.81 


68.78 


Average hospital length of stay (day) per patient 


4.87 


13.52 


8.96 


12.21 


10.51 
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(a) Patient-physician network 




(b) Corresponding PCN 



Figure 3 Construction of PCN from research dataset. The red circle represents physician and the gray triangle represents patient. 



its product with patient age. To show controlling effect, 
the product of network measure and patient age must 
show significant association with hospital outcome vari- 
ables in these regression models [55]. Out of these six 
models, this product shows a significant association in 
only two cases (i.e. the second and the third models of 
Table 4). That means patient age moderates only the 
relations of betweenness centralisation and density of 
the PCN with hospitalisation cost. In most cases (i.e. 4 
out of 6), patient age does not moderate the relation be- 
tween PCN attributes and hospital outcome measures. 
This can be explained by the fact that we do consider 
average age of all patients in calculating patient age for 
a PCN. On the other hand, studies of present healthcare 
literature consider patient age at the individual level, 
not at the aggregate level as like this study. 



Exponential random graph (ERG) model and physician 
collaboration network (PCN) 

Pnet a [1] has been used in this study to fit ERG models 
with different types of PCNs (i.e. low-cost versus high-cost 
and low-readmission versus high-readmission). After fol- 
lowing several iterative processes, the model (i.e. 2-star, 
3-star, alternating-k-stars, alternating- k-triangles and 
alternating-k-two-paths model) had been found to fit with 
PCNs. The results for this model are shown in Table 5. The 
weights of different micro -structures (e.g. 2-star and 3-star) 
of this model can be tested using t-value (also known as 
t-statistics), which is defined by dividing the estimate by its 
standard error. Thus, the t-value measures how many 
standard errors the estimate is away from zero. Generally, 
any t-value greater than +2 or less than -2 (i.e. absolute 
t-value is greater than 2) is acceptable. The higher the 



Table 2 Descriptive statistics (M indicates mean and STD indicates standard deviation) and pair wise correlation 
coefficient values for all variables (i.e. SNA variables, hospitalisation cost and readmission rate) used in this study 



[1] [2] [3] [4] [5] 



[1] Degree centralisation 










(M = 0.75 and STD = 0.13) 










[2] Betweenness centralisation 


-0.062 








(M = 0.25 and STD = 0.14) 










[3] Network density 


0.045 


-0.046 






(M = 0.27 and STD = 0.11) 










[4] Hospitalisation cost 


0.112 


-0.264* 


0.282" 




(M = 24009.9 and STD = 6783.3) 










[5] Readmission rate 


-0.373" 


-0.283" 


0.358* 


0.098 


(M = 11.64 and STD = 8.48) 











**. Correlation is significant at the 0.01 level (2-tailed). 
*. Correlation is significant at the 0.05 level (2-tailed). 
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Table 3 Linear regression models between each of network attributes (i.e. degree centralisation, betweenness 
centralisation and density) of PCN and hospital performance measures (i.e. hospitalisation cost and readmission rate) 



Model 


Dependent variable 


Independent variable 


R 2 value 


P 


Constant 


Significance 


1 


Hospitalisation cost 


Degree centralisation 


0.012 


5906.42 


19545.39 


0.309 


2 


Readmission rate 


Degree centralisation 


0.139 


-37.87 


45.48 


0.003 


3 


Hospitalisation cost 


Betweenness centralisation 


0.107 


-12384.79 


27101.96 


0.015 


4 


Readmission rate 


Betweenness centralisation 


0.112 


-25.18 


18.19 


0.010 


5 


Hospitalisation cost 


Density 


0.118 


17310.51 


19635.69 


0.009 


6 


Readmission rate 


Density 


0.196 


49.08 


3.427 


0.000 



t-value, the greater the confidence is shown by the param- 
eter under consideration as a predictor. Low t-value is the 
indication of a low reliability of the predictive power of 
that parameter [56]. To compare t-values of different 
networks, researchers utilise the t-test method [49]. 

The parameter interpretation of the fitted ERG model 
is summarised as follows. The positive 2-star parameter 
indicates that there is a tendency for multiple network 
partners. There is a significant difference in t-value of 
the 2-star parameter between high-cost PCNs and low- 
cost PCNs (Table 5). A t-test in Table 6 shows this 
significance (t (10) = 2.13, p<0.05). The results show 
that on average, the 2-star parameter for low-cost PCNs 
(M = 1.82, SE = 1.33) is less positive (in t-values) than 
the parameter for high-cost PCNs (M = 9.67, SE = 3.45). 
This indicates that the tendency for multiple network 
partners is more for high-cost PCNs than low-cost PCNs. 
It can be suggested from this trend that most of the 
actors of high-cost PCNs have multiple network partner- 
ships with others. This means that high-cost PCNs are 
well connected, but there is a low probability of having 
any network-hub (i.e. a highly connected actor), which 
indicates that these networks (i.e. high-cost PCNs) are 
decentralised. 

Table 4 Linear regression models for checking controlling 
network attributes (i.e. degree centralisation, betweenness 
measures (i.e. hospitalisation cost and readmission rate) 



Therefore, as PCNs become more centralised, the per- 
formance (i.e. inverse of cost) of the network will increase 
compared to decentralised PCNs. 

There is a significant difference in t-value of the 
alternating-k-star parameter between high- readmission 
PCNs and low-readmission PCNs (Table 5). A t-test in 
Table 6 shows this significance (t (10) = 1.75, p<0.05). The 
results show that on average, the k-star parameter for high- 
readmission PCNs (M= -3.69, SE = 0.87) is more negative 
than the parameter for low-readmission PCNs (M = -1.73, 
SE = 0.71). The negative alternating-k-star parameter indi- 
cates that networks with some higher degree nodes are 
less probable, which means there is no actor playing the 
role of network-hubs. This means that high-readmission 
PCNs are more decentralised. This indicates that as PCNs 
become more centralised, the performance (i.e. inverse of 
readmission rate) of the network will improve compared 
to decentralised PCNs. 

A significant difference in t-value for the alternating- 
k-2-paths parameter has been noticed between high- 
readmission PCNs and low-readmission PCNs (Table 5). 
A t-test in Table 6 shows this significance (t (10) =3.04, 
p<0.05). The results show that on average, the alternating- 
k-two-paths parameter for high-readmission PCNs (M = 

effect of patient age on the relation between each of 
centralisation and density) of PCN and hospital performance 



Model 


Dependent variable 


R 2 value 


Constant 


Independent variable 


P 


Significance 


1 


Hospitalisation cost 


0.102 


20016.04 


Degree centralisation 


-26621.75 


0.084 










Degree centralisation * Age 


463.16 


0.102 


2 


Hospitalisation cost 


0.227 


27106.44 


Betweenness centralisation 


-102698.68 


0.000 










Betweenness centralisation * Age 


1318.36 


0.000 


3 


Hospitalisation cost 


0.188 


19216.11 


Density 


-60258.64 


0.016 










Density * Age 


1155.44 


0.001 


4 


Readmission rate 


0.113 


45.82 


Degree centralisation 


-50.37 


0.078 










Degree centralisation * Age 


0.18 


0.622 


5 


Readmission rate 


0.094 


18.17 


Betweenness centralisation 


-75.92 


0.110 










Betweenness centralisation * Age 


0.745 


0.273 


6 


Readmission rate 


0.081 


1.318 


Density 


20.08 


0.677 










Density * Age 


0.314 


0.651 



Table 5 The results from high-order model (i.e. 2-star, 3-star, alternating-k-stars, alternating-k-triangles, alternating-k-two-paths model) 



Effects 


Estimates 


stderr 


Est./stderr 


Estimates 


stderr 


Est./stderr 


Estimates 


stderr Est./stderr 


Estimates 


stderr 


Est./stderr 


Estimates 


stderr 


Est./stderr 


















Low cost 




















N1 






N2 






N3 






N4 






N5 




2-star 


0.04 


0.01 


6.25 


0.12 


0.04 


2.66 


-0.04 


0.04 


-0.97 


0.03 


0.01 


2.02 


-0.02 


0.02 


-0.84 


3-star 


0.00 


0.00 


0.24 


0.00 


0.00 


0.86 


0.01 


0.00 


2.55 


0.00 


0.00 


-0.41 


0.00 


0.00 


3.03 


AS 


-1.95 


0.42 


-4.60 


-2.39 


1.24 


-1.92 


-1.24 


0.61 


-2.01 


-1.36 


0.38 


-3.58 


-3.55 


1.23 


-2.89 


AT 


1.28 


0.10 


12.64 


-0.05 


0.07 


-0.73 


1.37 


0.19 


7.04 


1.64 


0.14 


1 1.54 


1.43 


0.11 


12.68 


A2P 


-0.05 


0.01 


-5.70 


-0.17 


0.03 


-6.02 


-0.17 


0.02 


-7.98 


-0.02 


0.02 


-0.85 


-0.14 


0.01 


-15.40 


















High cost 




















N1 






N2 






N3 






N4 






N5 




2-star 


0.04 


0.00 


19.35 


0.03 


0.00 


13.71 


0.02 


0.00 


9.53 


-0.04 


0.03 


-1.34 


0.08 


0.01 


7.19 


3-star 


0.00 


0.00 


-3.84 


0.00 


0.00 


-1.74 


0.00 


0.00 


-1.10 


0.00 


0.00 


2.58 


0.00 


0.00 


-0.77 


AS 


-5.45 


1.03 


-5.28 


-5.05 


1.14 


-4.44 


-3.30 


0.55 


-5.96 


-9.60 


8.41 


-1.14 


-3.68 


1.30 


-2.83 


AT 


0.76 


0.06 


1 1.74 


1.06 


0.08 


12.72 


1.72 


0.08 


20.66 


-0.43 


0.06 


-6.59 


0.33 


0.09 


3.51 


A2P 


-0.05 


0.00 


-14.60 


-0.06 


0.00 


-13.60 


-0.04 


0.00 


-7.76 


-0.23 


0.02 


-11.99 


-0.07 


0.01 


-5.40 
















Low readmission rate 




















N1 






N2 






N3 






N4 






N5 




2-star 


-0.81 


1.78 


-0.45 


0.05 


0.03 


1.76 


0.01 


0.06 


0.21 


0.03 


0.02 


1.41 


0.02 


0.09 


0.16 


3-star 


0.08 


0.42 


0.18 


0.00 


0.00 


1.75 


0.01 


0.00 


1.89 


0.00 


0.00 


1.40 


0.01 


0.01 


1.06 


AS 


0.15 


2.98 


0.05 


-1.76 


0.77 


-2.28 


-5.05 


1.94 


-2.60 


-3.21 


0.89 


-3.61 


-0.19 


0.98 


-0.20 


AT 


1.31 


0.33 


3.93 


0.68 


0.12 


5.49 


0.70 


0.14 


4.94 


1.71 


0.19 


8.79 


-0.02 


0.11 


-0.20 


A2P 


0.17 


0.36 


0.46 


-0.1 1 


0.03 


-4.05 


-0.24 


0.02 


■13.12 


-0.1 1 


0.01 


-9.88 


-0.16 


0.06 


-2.76 
















High readmission rate 




















N1 






N2 






N3 






N4 






N5 




2-star 


0.02 


0.01 


0.05 


0.06 


0.01 


-0.03 


-0.02 


0.01 


0.02 


0.05 


0.03 


0.04 


0.30 


0.26 


-0.17 


3-star 


0.00 


0.00 


0.07 


0.00 


0.00 


0.05 


0.00 


0.00 


0.08 


0.00 


0.00 


0.04 


-0.02 


0.02 


-0.17 


AS 


-5.23 


1.29 


0.04 


-8.65 


1.25 


-0.02 


-6.65 


2.47 


0.05 


-7.54 


2.95 


0.01 


-26.01 


11.73 


-0.34 


AT 


0.48 


0.11 


-0.02 


0.97 


0.10 


0.01 


0.46 


0.12 


-0.01 


0.97 


0.15 


0.03 


-0.07 


0.07 


-0.34 


A2P 


-0.15 


0.01 


0.05 


-0.11 


0.00 


0.05 


-0.18 


0.01 


0.09 


-0.11 


0.02 


-0.04 


-0.33 


0.01 


0.09 
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Table 6 f-fesf for the t-values of different parameters of ERG model 

t-values of ERG model 

Parameter Low readmission High readmission Low cost High cost t-test P (T<=t) 

rate (mean) rate (mean) (mean) (mean) one-tail 

est/sd(2-star) 1.82 9.69 2.13 0.03 

est/std(AS) -1.73 -3.69 1.75 0.04 

est/std(A2P) -5.87 -21.58 3.04 0.01 



-21.58, SE = 4.55) is more negative than the alternating- 
k-two-paths for low -readmission PCNs (M = -5.87, 
SE = 2.47). The negative parameter of alter nating-k-two- 
path indicates that the network does not tend to form cy- 
cles and this tendency is higher in high-readmission PCNs. 
So in high-readmission PCNs, the tendency to form cycles 
will be less. 

Discussion and conclusion 

In this study, PCNs are constructed from the information 
of physicians' visits to patients during their hospitalisation 
period. It is assumed that collaboration emerges between 
two physicians when they visit a common patient. It is a 
standard professional practice that when physicians visit 
patients they give advice or suggestions to patients based 
on their health condition and previous medication history 
deposited in the patient log book. All previous advice or 
suggestions prescribed by any physician to a patient have 
been taken into consideration during any subsequent 
physician visit to that patient. In addition, physicians often 
have been informed about the patient condition by other 
physicians who previously visited that patient. This kind 
of practice culture in healthcare organisations or hospitals 
establishes the validity and reliability of the construction 
process of PCN, and the generic nature of the research 
findings. 

It is noticed that SNA measures of density for PCN has 
a positive correlation with hospitalisation cost and re- 
admission rate. In a dense PCN, an increased number of 
links exists among physicians. Although connections with 
peers enable physicians a faster sharing of known know- 
ledge [57], links with many peers significantly impacts an 
individuals opportunity to share or create knowledge in a 
network (e.g. PCN). This is because when an individual 
has many links in a network, she needs to spend more 
time to maintain these relationships. Moreover, she will 
receive repetitive or contradictory knowledge, in addition 
to new knowledge, from many individuals which will make 
it difficult to summarise these shared knowledge [58]. Un- 
like density, the SNA measure of betweenness centralisa- 
tion has negative correlation with both hospitalisation cost 
and readmission rate. From the perspective of a PCN 
structure, a high betweenness centralisation indicates that 
the structure of the corresponding PCN follows a star-like 
or centralised structure since betweenness centralisation 



reaches its highest value of 1 for a star network. A star- 
like or centralised network has few actors with higher 
betweenness centrality values. In this type of network, only 
a small number of actors play major collaboration and 
communication role. Therefore, in their corresponding 
hospitals, healthcare managers or administrators have to 
encourage or establish a star-like or centralised PCN in 
order to reduce both hospitalisation cost and readmission 
rate. A PCN with a flat network structure (i.e. members of 
that PCN have almost equal network participation) will 
have high hospitalisation cost and readmission rate. 

Although this study finds that social network measures 
have statistically significant correlations with hospitalisa- 
tion cost and readmission rate the corresponding correl- 
ation coefficient values do not show perfect correlations 
(i.e. a correlation coefficient value of 1) among them. The 
correlation coefficient values are ranging from 0.112 to 
0.358 in absolute values. That means none of the relations 
shows perfect correlation. However, five of these values 
(see Table 2) are statistically significant at 0.01 and 0.05 
levels (2-tailed). This is because of the sample size used in 
this study. This study uses 85 PCNs to explore the effects 
of different SNA measures on hospitalisation cost and re- 
admission rate. A small correlation coefficient value could 
be statistically significant if sample size is high; whereas, 
for a small sample size a high correlation coefficient value 
would not be statistically significant [59]. A correlation co- 
efficient value of 0.04, for instance, would be statistically 
significant for a sample size of 10,000 [59]. We also quan- 
tify the impact of social network measures on hospital 
outcome variables using simple linear regression models 
(see Table 3). 

From the ERG model (i.e. 2-star, 3-star, alternating- 
k-stars, alternating- k-triangles and alter nating-k-two-paths 
model), significant differences are noticed in t-values for 
different micro- structures between high-readmission PCNs 
and low -readmission PCNs. The alter nating-k- star param- 
eter shows a stronger negative value for high-readmission 
PCNs compared to low -readmission PCNs. Similarly, 
alternating-k-two-path shows more negative values for 
high-readmission PCNs compared to low-readmission 
PCNs. In summary, high-readmission PCNs are attributed 
with tronger negative values for alternating-k-star and 
alternating-k-two-path parameters. Negative t-value for an 
alternating-k-star parameter implies that networks with 
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well-connected nodes are less probable (i.e. lack of 
the presence of network-hub). Negative t-value for an 
alternating-k-two-path parameter reveals that network 
actors are less likely to form cycles (i.e. networks are most 
sparse). From the interpretation of the findings of these 
two parameters (i.e. alternating-k-star and alternating-k- 
two-path) it can be concluded that high-readmission PCNs 
are decentralised. And low-readmission PCNs are more 
centralised compared to high-readmission PCNs. Low-cost 
PCNs and high-cost PCNs also have micro- structural dif- 
ferences that are statically significant. The t-value of 2-star 
parameter is more positive for high-cost PCNs compared 
to low-cost PCNs. So, in high-cost PCNs, actors tend to 
have multiple partnerships with other network actors. 
That means high-cost PCNs are well connected and have a 
low chance of having any network-hub (i.e. highest degree 
actors). Therefore, high-cost PCNs are less centralised. 
And low-cost PCNs are more centralised compared to 
high-cost PCNs. To summarise, the findings from ERG 
models, both low-readmission rate and low-cost PCNs are 
more centralised compared to their counterparts. There- 
fore, in their corresponding PCNs, physicians have to be 
close to each other; and they should not work standalone 
or in disconnected small groups. A centralised PCN en- 
ables effective knowledge sharing among its member phy- 
sicians, which eventually leads to better patient care [57] . 

In respect of the network data analysis using ERG 
model, this study intends to explore micro-level struc- 
tures (e.g. 2-star and 3-star) that are associated with dif- 
ferent PCNs characterised by the highest hospitalisation 
cost versus lowest hospitalisation cost and highest re- 
admission rate versus lowest readmission rate. For this 
purpose, this study considers only the top 5 PCNs from 
all these four groups instead of considering all 85 PCNs 
of our research dataset. 

This research is not without its limitations. First, we 
test relations between PCN attributes and healthcare 
care outcome measures and ERG models using the 
health insurance dataset only for THR patients. Thus, 
we need to consider dataset for other patients such as 
knee surgery patients or patients suffering from brain 
cancer in order to claim the general nature of the find- 
ings of this study. Second, we consider only quantitative 
measures (i.e. hospitalisation cost and readmission rate) 
as outcome variables. We do not consider any qualitative 
measures (e.g. patients' satisfaction) as outcome variables. 
Finally, we consider only 20 PCNs for ERG modelling, 
which significantly limits the interpretation of the ERG 
findings. This is because, in this study we aim to explore 
structural differences in PCNs classified as high and low 
in terms of hospitalisation cost and readmission rate. 

To conclude, this study first proposes a way to capture 
networks that evolve among physicians during the course 
of providing treatments to hospitalised patients. Second, 



SNA measures are utilised to explore PCNs. It is noticed 
that density has positive correlation with hospitalisation 
cost and readmission rate; whereas, betweenness central- 
isation is negatively associated with hospitalisation cost 
and readmission rate. Degree centralisation shows no sig- 
nificant correlation with hospitalisation cost and negative 
correlation with readmission rate. Finally, an ERG model 
is fitted with different types of PCNs (i.e. low-cost versus 
high-cost and low-readmission versus high- readmission). 
From the ERG model, it is found that PCNs, which are at- 
tributed with less negative t-values for alternating-k-star 
and alternating-k-two-path parameters, and lower t-values 
for 2-star parameter, are more conducive to performance 
in terms of low hospitalisation cost and low readmission 
rate for patient hospital admissions. Healthcare managers 
and hospital administrators may follow the findings of this 
study in promoting the physician collaborations structure 
within their organisations. 

Endnote 

a http://www.sna.unimelb.edu.au/pnet/pnet.html. 
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